Triplet pairing and upper critical field in the mixed state of d-wave superconductors 
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We show that an additional triplet component of the order parameter is generated in the vortex 
phase of the d-wave superconductor. Spatial variations of the triplet component are analyzed for 
a strong spin-orbit coupling. Corrections to the London equation and an unusual temperature 
dependence of the upper critical field, H C 2(T), are obtained in the case of the weak spin-orbit 
coupling. 
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There is a common belief that the copper-based high- 
temperature superconductors have nontrivial order pa- 
rameter (OP) transforming as B\ g {x 2 — y 2 ) representa- 
tion of Dih point group [l[ . In the framework of the BCS 
theory this " d-wave" symmetry is related to the 'internal' 
coordinate of the OP, r = ri — r 2 , assuming its decompo- 
sition as A(ri,r 2 ) = A (R) A (r), where R = (ri+r 2 )/2 is 
the center-of-mass coordinate This symmetry of the 
OP implies that the spectrum of elementary excitations 
is gapless in certain directions so that the d-wave gap is 
very sensitive to an external perturbation. For example, 
applying a uniform magnetic field generates a secondary 
singlet component of the OP, that lowers the symmetry 
and opens a full gap in the nodal directions [!, 0] . P-wave 
component of the OP could be also generated by a surface 
induced spin-orbit coupling (SOC) ||. Similarly, super- 
conducting currents give rise to ip secondary component 
and partial opening of the gap in the nodal directions 
@. Appearance of ip component 0] is associated with 
the Lifshitz invariants (LI) in the free energy (FE) 
More recently it has been shown that the triplet OP iA t 
is generated in the mixed state of type-II superconduc- 
tors @ . It was argued that this effect is a consequence of 
broken symmetry in the spin space due to the paramag- 
netic effect and a broken translational invariance. This 
effect was called as type-IV supeconductivity Q. 

Phcnomcnologically nontrivial secondary OP is gen- 
erated by linear coupling of the primary OP gradients 
to the secondary OP. From the microscopic point of 
view the interaction of electrons on the Fermi surface 
is repulsive in the channel of secondary OP. Neverthe- 
less, non-diagonal terms corresponding to the coupling 
to the primary OP leads to nonzero value of the sec- 
ondary OP[6|, |8j. The relative amplitude of the p-wave 
component depends on the strength of SOC and may be 
of the order 10-20% @. 

Previously it was shown that the superconducting cur- 
rent in the d-wave superconductor reduces the symmetry 
group and generates the secondary triplet OP @. It is 
possible since the symmetry allows for the presence of 
terms in the Ginzburg-Landau FE functional which are 
of the first order in gradients. Usually these terms are 
allowed for crystals where the inversion symmetry is bro- 
ken. Inversion is the symmetry operation of group 
and generation of the LI is possible if the secondary com- 



ponent breaks the inversion symmetry 0, Q . Here we in- 
vestigate the electromagnetic response in the mixed state 
of the type-II d-wave superconductor in the presence of 
LI @j 0] and formulate the criteria for experimental ob- 
servation of the effect. 

First let us briefly discuss the symmetry properties of 
p-wave superconducting OP [10(. The triplet OP has 
9 components. In the case of the strong SOC spin is 
coupled to the lattice and transforms together with the 
lattice. P-wave component of D^h point group is written 
as d(k, R) = (p x (R)k x + p y (IV)k y )z,, where z is a unite 
axial vector in the spin space and ^(R) and p y (R) plays 
the role of the OP in the Ginzburg-Landau functional. 
In the case of the weak SOC the spin and the orbital 
parts of the OP transforms independently and d(k, R) = 
p x (R)k x + Pj,(R)/cj,, where two vectors in the spin space 
p x (R) , p y (R) play the role of the OP. In the latter case 
the coupling between singlet and triplet OP is possible 
only in the presence of the magnetic field. 

Strong SOC. The FE of d-wave superconductor per 
unite of volume can be written as: 
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where ip and p x%v are the primary d-wave and the sec- 
ondary p-wave OP, D = — iV — 2eA, A is the vector po- 
tential for magnetic field H = V x A , a, a p , (3 > 0, n are 
real constants. Here we write explicitly the LI describ- 
ing p-wave component. In Eq.(l) the additional term, 
quadratic in p x , p y ,is written to secure p x — p y — solu- 
tion in the uniform state ipd = const. 

Minimizing the FE for variation of Sip, 5p x , 8p y we ob- 
tain the following set of equations: 
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Substituting Eq.(3) to Eq.(2) we obtain standard 
Ginzburg-Landau equation 
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: DV + a(T - T c )ip + f3\ip\ 2 ip = 
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with renormalized effective mass (m*) -1 = mT x — 
2r] 2 /a p . The renormalization of the effective mass was 
derived in Eq.(7) of Ref.[6]. If renormalized effective 
mass is negative (to*) -1 < (2rj 2 /a p > to -1 ) the heli- 
cal phase is formed [7[ and higher order gradient terms for 
secondary OP should be included to the FE, Eq.(l). The 
properties and the thermodynamics of helical phases have 
been studied in Ref.Q and are not considered here. If 
renormalized effective mass (to*)" 1 > {2rf/a p < to" 1 ) 
the FE, Eq.(l), describes the formation of the secondary 
OP in the presence of the magnetic field or in the presence 
of the supercurrent [6j . Therefore we assume in the fol- 
lowing that (to*) -1 > and higher order gradient terms 
are not essential in generation of the secondary OP. 

Minimizing the FE for variations of the vector poten- 
tial SA and substituting Eqs.(3,4) we obtain supercon- 
ducting current: 
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Now let us calculate the OP and the properties of 
the mixed state in the external magnetic field. It is 
convenient to introduce dimensionless variables, i/> = 
■0o/, Px, v = ipoVx.y, where Vo = VH T = T c)\/P, 



X(T)p, X(TY 
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£(T)A, where 



£,{T) 2 = - 2m , a ( T _ Tc) and h = |§£(T)A(T)H. Separat- 
ing the phase and the modulus of the OP, f = fo exp (<fi) , 
and redefining the vector potential the system of equa- 
tions is reduced to the equation for the modulus of the 
OP f — fo and for the magnetic field h[ll[: 
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where n = A/£ is Ginburg-Landau parameter. Amplitude 
of the p-wave component is expressed in terms of fo and 
h as: 
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Therefore the problem is reduced to the standard solution 
of the Eqs.(6,7) and then to the calculation of the p-wave 
amplitudes applying Eq.(8). 

Weak magnetic field H c \ < H <C H c i- Single vortex. 
Solution of Eqs. (6,7) are known and can be written in 
the form fill: 



fo = cp + 
h = h(0) - cp 2 /2n + 
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where c is the constant of the order of 1 . This expansion 
is valid near the vortex core p <^ 1/k. The amplitude of 
the p-wave component in that case is: 



•qc 



ly 



(10) 



This effect is consistent with the symmetry arguments 
and was briefly discussed in connection with superfluidity 
of 3 ife[I3]. Near the vortex core the screening current is 
flowing around the vortex as well as the gradient of the 
modulus of the OP is finite. As it was shown in Rcf. 6] 
currents generate ip secondary OP. On the other hand a 
real vector breaks the inversion symmetry and generates 
the real p-component of the secondary OP [5] . Far from 
the vortex core but in the region where the screening 
current is strong 1/k <C p <C 1 fo~ 2 dh/dp = — (up)" 1 

and fo = 1 - ( K PY 
component are 
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This is again consistent with the previous results. Far 
from the vortex core the modulus of the OP is almost 
constant (gradient is small) and the screening current is 
strong. Therefore only ip secondary OP survives in that 
area[6j . 

Strong magnetic field (H C 2 — H) / H C 2 < 1 Magnetic 
field in that case is determined by the formula ll( : 
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where e = K — ZJ Hc2 . Spacial dependence of f is deter- 



mined by[TT| : 
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where |c | 2 = (H - H c2 )/f3H c2 and /3 = 1.1596. Substi- 
tuting Eqs. (12,13) to Eq.(8) we obtain 

^ = -^ = ^(§+#)- d4) 
a p !;K ox ay 

Weak SOC. In the strong magnetic field when Zee- 
man energy is larger then the energy of the SOC spin 
decouples from the lattice. In that case vector func- 
tions p x (R) , p y (R) will be coupled to the magnetic field 
H(similar effect was considered in Ref . [l3| ) . Therefore 
linear in gradients term in the FE has the following form: 
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Here we use cartesian basis in the spin space [lij and 
therefore the triplet OP has zero projection of the spin 
to the direction of the magnetic field [§| . Minimizing the 
FE for variation of the OP Sip, Sp x , Sp y we obtain the 
following equation for ip: 

^T) 2 4,+a{T-T c )^+m 2 i>-— HD 2 (W) = 0, (15) 

2m a.p 

where p x>y — ±^-D x , y (H.ip). If we assume that the mag- 
netic field is constant and equal to the external magnetic 
field Eq. (15) is reduced to Eq. (4) with the field de- 
pendent effective mass (m*(i/)) _1 = to" 1 — 2rf 2 H 2 /a p . 
Straightforward variation with respect to vector poten- 
tial A gives the expression for superconducting current 
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To simplify this equation let us assume that the magnetic 
field is a constant. After simple calculations we get the 
formula for j s , 
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where v = -^j^(a(T ~ T c )|^| 2 + f |V| 4 ). It is seen 
from this equation that the correction to the supercur- 
rent is similar to the ordinary Hall current, where an 
effective electric field is proportional to the gradient of 
the superfluid density, n s = |i/>| 2 . 

When the magnetic field increases further effective 
mass m*(H) becomes small and next order terms in the 
gradients should be considered in the expansion of the 
FE density. This is important for calculation of the 
upper critical field H c i- For the calculation of H C 2 we 
analyze linearized equation with constant magnetic field 
directed along z axis. Therefore only z component of 
p x ,y OP is relevant. We also write second order gradi- 
ent terms for p-wave component in the FE density as 



|Dpj,| 2 ). Introducing new dimensionlcss 



variables, H = 27r |° H , £ 2 = 2 maT ' ®° 1S * ne nux q uan ta 
and x — > x\^H/^q, one obtains the following set of equa- 
tions: 

H(-ij}" + x 2 ^) +t^ + H 3/2 v(p' x + ixp y ) = 
H{-p" x + x 2 Px ) + bp x - vcV? 12 ^ = (18) 
H(-p'y + x 2 p y ) + b Py - ivcH s/2 xip = 0, 



where v — 2 ^aT ' c = to p /to ~ 1 and b — a p c/aT c ~ 
1, t = (T — T c )/T c . Here we use Landau gauge A = 
(0, Hx, 0) and therefore assume that all functions depend 
on x only. Using standard substitution: 

ip = s n cxp (-X 2 /2)H n (x) 
P*,y = Y,1n' y) ^p(-x 2 /2)H n (x), (19) 

where H n {x) are hermitian polynomials. After simple 
calculations we obtain the following algebraic equations 
which determines critical temperature in the magnetic 
field: 
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Upper critical field is determined by the highest possible 
critical temperature which corresponds to n — in the 
Eq. (20). Hence the critical transition temperature in 
the magnetic field is determined as 



T C (H)/T C = 1-H 



v 2 cV? 
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There are a number of interesting observations based 
on Eq.(21). The presence of LI increases critical tem- 
perature of the superconductor. This effect is well 
known Q. More interestingly, the correction to the up- 
per critical field becomes nonlinear and leads to the pos- 
itive curvature of H c2 as the function of temperature: 
d 2 H C 2(T)/dT 2 > 0. Further increase of the magnetic 
field may lead to the recovery of superconductivity in 
higher fields [TB|. It is also consistent with the previous 
results Q, where it was shown that increase of the criti- 
cal temperature due to formation of inhomogeneous he- 
lical phase is proportional to rj 2 . Therefore, calculated 
upper critical field should be substantially larger then 
paramagnetic or Clogston limit [H. flrl| . Generalizing this 
result to the weak SOC we can see that the shift of the 
critical temperature due to purely paramagnetic effect 
is proportional to j] l2 H 2 . This is exactly the results for 
T C (H) which follows from Eq.(21). We have to remember 
that in high magnetic field higher order gradient terms 
should be considered. The higher order orbital effect is 
proportional to Ti 2 and reduce the effect of LI. If higher 
order orbital effects are strong enough the reentrant su- 
perconductivity in high magnetic field will be suppressed 
completely. 

In Figure 1 we plot the temperature dependence of 
H C 2 . To avoid infinite increase of critical temperature in 
the high field we take into account higher order gradient 
terms in the formula Eq.(21). As a result critical field is 
determined by the formula: 
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FIG. 1: Temperature dependence of the upper critical field 
H C 2, parameter v 2 c =1, 1.1, 1.12, 1.13, 1.15 

where 7, j p ~ 1 dimensionless constants describing the 
effect of the higher order gradient terms in the Landau 
expansion. To demonstrate the temperature dependence 
of the critical field we choose 7 = 1/4, b — 1, -f p = 1 and 
v 2 c is changing from 1, to 1.15. The results, presented in 
Fig. 1 are restricted by standard range of applicability 
of Ginzburg-Landau theory and the limit T — > should 
be considered only as a qualitative. In that case more 
accurate microscopic equation [13] should be analyzed. 

To estimate the strength of the effect we refer to 
Gorkov's equation formulated for the case of the strong 
SOC inpl], and in the case of paramegnetic effect injal. 
It follows from the analysis in the strong SOC limit [6[ 
A p /A d m Aso/Ef w A S o^bH c2 / A 2 d , where A Pjd is the 
p(d) gap, respectively, A so is the energy of the SOC, Ej 
is the Fermi energy, and is the Bohr magneton[l8(. 
From the phenomenological consideration it follows that 
Ap/Ad ~ v/^p- Therefore LI are determined by the 



SOC r]/£a p w Aso/Ef. In the case of the weak SOC 
the estimate reads A p /Ad w r}'H/£a p « hbH&Z Ad\^. 
In that case LI are determined by by the paramagnetic ef- 
fect. Therefore at the value of the field H ss iJ c2 Aso I Ad 
we expect a crossover from the strong to the weak SOC 
limit. The characteristic value of triplet gap was evalu- 
ated in[8( and may be of the order of 0.8 of in super- 
conductors where H C 2 is close to the paramagnetic limit. 
As a result we can estimate the value of the parameter 
v 2 c w (H C 2/Ad) 2 ~ 1 in Fig. 1. Therefore we believe 
that the effect may be observed in materials with weak 
SO coupling and with high upper critical field if C 2> that 
exceeds paramagnetic limit. Possible candidates for this 
effect are organic superconductors, where highly nonlin- 
ear temperature dependence of the upper critical field is 
observed (l9| . On the other hand, the p-wave component 
of the order parameter in the mixed state of supercon- 
ductors may be observed by the direct vortex imaging 
with the scanning tunnelling spectroscopy. This effect 
should be strong in materials with the strong SOC. 

In conclusion we have shown that the additional triplet 
secondary OP is generated in the vortex (Abrikosov) 
phase of type II d-wave superconductor. This effect is 
caused by the superconducting currents and by the spa- 
tial variations of the OP in the presence of the magnetic 
field. According to estimates using the microscopic the- 
ory the triplet component may be as large as 80% of the 
primary d-wave component in the field near H C 2 @ . Non- 
linear corrections in the vector potential to the London 
equation are derived. The upper critical field H C 2{T) is 
calculated in the limit of weak SOC. It is shown that tem- 
perature dependence of H C 2(T) is highly nonlinear near 
the transition to the helical phase. 
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